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Black hole perturbation theory is typically studied on time surfaces that extend between the 
bifurcation sphere and spatial infinity. From a physical point of view, however, it may be favorable 
to employ time surfaces that extend between the future event horizon and future null infinity. This 
framework resolves problems regarding the representation of quasinormal mode eigenfunctions and 
the construction of short-ranged potentials for the perturbation equations in frequency domain. 



I. INTRODUCTION 

Black hole spacetimes are typically represented in co- 
ordinates in which the time hypersurfaces intersect at the 
bifurcation sphere and at spatial infinity. Because of the 
simplicity of these coordinates, black hole perturbation 
theory is typically studied along such time surfaces [T}]1] . 

In this report, I argue that it is favorable — especially 
for numerical computations — to employ time surfaces 
that extend between the future event horizon and fu- 
ture null infinity. These surfaces avoid the pathological 
behavior at the bifurcation sphere and at spatial infinity. 
They are called (future) horizon-penetrating and (future) 
hyperboloidal (Fig. [I]). 

The physical argument is fairly concise. Consider the 
standard representation of the Schwarzschild metric 



singularity 




g = -fdt^ + f-Ur' 



r2 dn^. 



with / := 1 — 



2M 



Here, r is the areal radius, dfl^ is the standard metric 
on the unit sphere, and M denotes the mass of the black 
hole. The surfaces of constant time t intersect at the 
bifurcation sphere, B, and at spatial infinity, i'^ (dashed 
lines in Fig.[T]). However, an astrophysical black hole that 
forms by gravitational collapse does not possess a bifur- 
cation sphere, and observers of gravitational radiation 
do not have access to spatial infinity. From this point of 
view, it is appealing to study black hole perturbations on 
spacelike hypersurfaces extending from the future event 
horizon, H^, to future null infinity, ^7+ (solid lines in 
Fig.[l]). 

Such hypersurfaces can be constructed using the height 
function technique [5 . Introduce a new time coordinate 



T{t,x')=t-h{x'), 



(1) 



where the height function, h, depends on spatial coor- 
dinates only. This transformation leaves the timelike 
Killing field in the exterior domain invariant. Observers 
at rest with respect to the black hole have the same rep- 
resentation in the new time coordinate. Many useful 
spacelike foliations of Schwarzschild spacetime, such as 
Eddington-Finkelstein or Painleve-Gullstrand foliations, 
have the above form. 

It is convenient to employ the tortoise coordinate 
= r + 2Mlog[r/(2M) — 1]. The event horizon in this 
coordinate is located at —oo. We choose a height 



FIG. 1: Penrose diagram of the domain of outer com- 
munications in Schwarzschild spacetime depicting standard 
Schwarzschild time surfaces (dashed lines) and future horizon- 
penetrating, future hyperboloidal time surfaces (solid lines). 



function that depends only on and define H := dh/dr^, , 
referred to as the boost function. The Schwarzschild met- 
ric becomes 

g = f {-dr'^ - 2HdTdr^ + (l - H"^) drl) + r^dfl^. 

The boost function for horizon-penetrating, hyper- 
boloidal time surfaces satisfies 



\H\ < 1, 



hm H 



±1. 



dH 
hm - — 

,->-±oo dr^ 



0. 



(2) 



These conditions imply that, asymptotically, the boost 
of r-hypersurfaces with respect to i-hypersurfaces ap- 
proaches the speed of light. An example is the boost 
function of constant mean curvature surfaces [51 [7] 



H = 



J 



J := 



Kr 



c_ 



(3) 



with a positive mean curvature, -fC > 0, and C > 
8M^K/3 (solid lines in Fig. [l] with C = 10 and X = 2). 

This physically motivated geometric framework has al- 
ready been applied in numerical studies of gravitational 
black hole perturbations in the time domain [8 10 . In 
the following, I present a theoretical discussion in the 
frequency domain. 
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II. REGGE-WHEELER-ZERILLI EQUATION 



The asymptotic form of this equation reads by ^ 



The radial Regge- Wheeler- Zerilh (RWZ) equation de- 
scribes metric perturbations of Schwarzschild spacetime 
via multipoles of a Fourier-decomposed master function 
[TJUj. Dropping the multipole indices, the radial RWZ 
equation can be written as 



dr'i 



[/ * 0. 



(4) 



The potential U depends on whether the perturbations 
have even or odd parity. Without loss of generality, we 
consider odd parity perturbations 

where ^ is the angular momentum number. The potential 
vanishes both at the horizon and at infinity. The asymp- 
totic behavior of the master function in these domains 
is 



C2e 



as 



±oo. 



(5) 



We arc interested in solutions that are purely ingoing 
near the horizon and purely outgoing near infinity, which 
translates into ^I* g-"^'". near the bifurcation sphere 
and ^' — >■ e'"^' near spatial infinity. Solutions with this 
asymptotic behavior describe exponentially damped os- 
cillations of the black hole, called quasinormal modes 
(QNMs) [IIHni- 

The frequency ui has, in general, both a real and an 
imaginary part. As a consequence, QNM eigcnfunctions 
blow up exponentially near the black hole and near infin- 
ity. This behavior is due to the pathological properties 
of the bifurcation sphere and of spatial infinity: The rep- 
resentation of the physical boundary conditions in these 
domains is unphysical. 

Switching to horizon-penetrating, hyperboloidal hy- 
persurfaces resolves this problem. The time transforma- 
tion ([T]) for the radial RWZ equation in frequency domain 
is equivalent to the rescaling of the master function 



(6) 



This rescaling, with a specific choice for h, is used by 
Dolan and Ottewill in their expansion method, which has 
been applied to the calculation of QNMs and Regge poles 
in the eikonal limit in Schwarzschild and Kerr spacetimes 
[131 [H] ^ ■ It is remarkable that they construct the trans- 
formation based on a very different reasoning. 

The transformation Q leads to the following equation 



(f d 

— + 2iuH— + Lu'^n - H^) + iuH' -U]ip^O. 
dri dr* 



(7) 



(8) 



with solutions 



i/) Cie 



as 



±CXD. 



The physical boundary condition with respect to the new 
foliation is simply that the solution is of order unity both 
at the horizon and at infinity, that is, ip — at both 
asymptotic ends. As a consequence, the rescaled eigcn- 
functions are regular in the infinite domain extending 
from the event horizon to null infinity. This feature ex- 
tends in a straightforward way to perturbations with even 
parity and arbitrary spin. 



III. TEUKOLSKY EQUATION 

The radial Teukolsky equation describes curvature per- 
turbations of Kerr spacetimes |4j . A practical difficulty in 
numerical computations with the radial Teukolsky equa- 
tion is its long-ranged potential. Among the various ef- 
forts to construct an equivalent equation with a short- 
ranged potential |16H19) . the most common approach in 
numerical calculations relies on a generalization of the 
Chandrasekhar transformation [20] devised by Sasaki and 
Nakamura [HI 112] • 

Sasaki and Nakamura reformulated the Teukolsky 
equation in the RWZ form, first for Schwarzschild space- 
time ^1] and then for Kerr spacetime [2^. Their for- 
malism was restricted to the gravitational case with 
spin s = —2. Hughes later generalized the formalism 
to arbitrary spins in Kerr spacetime [23'. The Sasaki- 
Nakamura transformation has been applied in various 
numerical calculations concerning the motion of particles 
in Schwarzschild and Kerr spacetimes [^IMT] . 

The framework presented in the introduction leads to a 
short-ranged potential with a simple transformation that 
has a clear physical and geometric interpretation. We 
discuss its application to the Teukolsky equation first in 
Schwarzschild spacetime, and then in Kerr spacetime, for 
arbitrary spins. 



A. Schwarzschild spacetime 

The radial Bardeen-Press-Teukolsky (BPT) equation 
in Schwarzschild spacetime for perturbations of spin 
weig ht s reads [3 14] 



d^ 2{{r-M){l + s)-M) d 2 



+ 



^ I thank Emanuelc Bcrti for pointing out these references. 



with the potential 
2iu)s 



-t- 0, 
(9) 



U ■ 



I 



(r - 3M) + + 1) - s{s + 1)). 
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The problematic term in Q is the first term of the po- 
tential that does not vanish at the horizon and falls off 
as r^^ near infinity. To construct a short-ranged po- 
tential, one that vanishes at the horizon and falls off at 
least as we rescale the master function according 
to its asymptotic behavior. The perturbations fall off as 
^-(2s-i-i) jjgj^j. infinity, and as near the black hole [1]. 
A further rescaling transforms the equation to a horizon- 
penetrating, hyperboloidal foliation as in ^. We set 



_ ^-(2s+l) 



r 



(10) 



The modified BPT equation reads 
2s 



:{r- M) + 2iujH 



dAp_ 



- H^) + iuH' - J7) V = 



(11) 



with the potential 
^ 2isuj 



{Tf{l-H)-M{l + H)) 



f / 2 M 

-f4 {e{i + i)-sis + i) + — (s + i) 

V r 



.(12) 



The asymptotic behavior of the boost function given in 
([2]) implies that this potential falls off as / near the hori- 
zon and as r^^ near null infinity. Therefore, the modi- 
fied potential is short-ranged. Note that all lowest order 
terms in (11) vanish at the asymptotic ends. 



B. Kerr spacetime 

The radial Teukolsky equation in Boyer-Lindquist co- 
ordinates reads 



dr \ dr 



where A := r'^ + — 2Mr, and 



U 



if 2 _ 2is{r - M)K 
A 



[/* = 0, 



+ AisLor — A, 



with K :— (r^ + a^)uj — ma, and 

A := Sim ~ 2amuj + a^w^ — s{s + 1). 

Here, Sim is the eigenvalue of the spherical harmonic [53] 
(it becomes Sim — ^(^ + 1) m the Schwarzschild limit). 

We write the Teukolsky equation in the tortoise coor- 
dinate to make the connection to the previous sections. 
The tortoise coordinate in Kerr spacetime is defined via 



A 



■ dr. 



The Teukolsky equation in the tortoise coordinate reads 
(r2 + a2)2^2^ 2G A_ 



dr? r^ dr* 



rfJ* = 0, 



where 

G = M{a^ - r^) + (r - M){r'^ + a^){s + 1). 

A geometric transformation that leads to a short- 
ranged potential can be given as 



(13) 



This transformation is very similar to ( 10 1 with one es- 



sential difference: We transform also the angular coordi- 
nate (j). The Boyer-Lindquist angular coordinate (p leads 
to pathologies near the horizon [55] . To cure this behav- 
ior, we introduce the angular coordinate </> defined by 



rdr* 



(14) 



This azimuthal transformation is used already by Teukol- 
sky in his seminal paper on rotating black hole perturba- 
tions to derive the boundary conditions near the horizon 
[3]. It is commonly used in numerical computations in 
time domain codes [551 1^ . 

The Teukolsky equation for the transformed eigenfunc- 
tion tp from (|13|) reads 



(r^ + a2)2 (p^ 2G dip C/ , 

with 

G = a^A+ir'^+a^) [rs{r - M) - ir{{r'^ + a^)u}H + ma)) , 
and 

U = 2isujr^ {rA{l- H)- M{r^ -a^){l + H)) 

—2iarA[m + aujH) + A (2a^ — A - 2Afr(s + 1)) 
-2maujr'^{r'^ + a^){l + H) 

The lowest order term U/r^ vanishes as A at the hori- 
zon and falls off as r^^ near null infinity due to the 
asymptotic behavior of the boost function ([2]). So the 
rescaled Teukolsky equation in horizon-penetrating, hy- 
perboloidal coordinates with the angular coordinate p of 
( 14 1 has a short-ranged potential. 



IV. DISCUSSION 

The framework described in this report has advantages 
over standard methods in black hole perturbation theory. 
The QNM eigenfunctions have a finite representation in 
an infinite domain, and the radial Teukolsky equation be- 
comes short-ranged with the transformation (13), which 



is simpler than the Sasaki-Nakamura transformation, and 
has a clear physical and geometric interpretation. 
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This method may lead to efRcient numerical codes in 
the frequency domain, especially in combination with 
compactification of the radial coordinate. Compactifi- 
cation would remove the cumbersome asymptotic expan- 
sions and the studies of inner and outer boundary effects 
in the numerical integration of ordinary differential equa- 
tions appearing in frequency domain computations. This 
technique should be studied in future research on astro- 
physically motivated examples, such as calculations of 
extreme mass ratio inspirals using the Teukolsky formal- 
ism, or the self-force approach. 

The geometric origin of this framework is a clear ad- 
vantage. It implies that one does not need to apply 
the transformation ( 13 1 back and forth in solving the 



Teukolsky equation numerically, as it is done with the 
Sasaki-Nakamura transformation . One can di- 

rectly solve the transformed equation, and extract the 
physical quantities, such as energy fluxes, at null infinity 
from the rescaled master function without going back to 
the Teukolsky equation in Boyer-Lindquist coordinates. 



The new framework can also be useful in analytic stud- 
ies. Dolan's expansion method as applied to Kerr space- 
time is restricted to equatorial and polar modes [TS] . The 
geometric framework presented in this report indicates 
that the angular transformation ( 14 1 may lead to an ex- 



tension of the expansion method to general modes in Kerr 
spacetime. It would also be interesting to study how this 
framework relates to Green function based regularization 
techniques of divergent integrals for the Teukolsky equa- 
tion with sources l30H32l. 
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